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The quantum-dynamical mechanism of photoinduced subpicosecond exciton dissociation and the 
concomitant formation of a charge-separated state at a TFB:F8BT polymer heterojunction is elu- 
cidated. The analysis is based upon a two-state vibronic coupling Hamiltonian including an ex- 
plicit 24-mode representation of a phonon bath comprising high-frequency (C=C stretch) and low- 
frequency (torsional) modes. The initial relaxation behavior is characterized by coherent oscillations, 
along with the decay through an extended nonadiabatic coupling region. This region is located in 
the vicinity of a conical intersection hypersurface. A central ingredient of the analysis is a novel 
effective mode representation, which highlights the role of the low-frequency modes in the nona- 
diabatic dynamics. Quantum dynamical simulations were carried out using the multiconfiguration 
time- dependent Hartree (MCTDH) method. 



The photophysics of 7r-conjugated organic semiconduc- 
tor systems is a key ingredient in the technological de- 
velopment of optoelectronic devices. Due to their spa- 
tially extended structure, these systems exhibit both the 
molecular characteristics of their components and the col- 
lective electronic excitations (exciton states) characteris- 
tic of lattice structures. The coupling between electronic 
and nuclear motions results in the localization ("self- 
trapping" ) of exciton states, along with the coherent nu- 
clear dynamics and vibronic coupling phenomena known 
from molecular systems. 0, H, UllEl H 0, H Experiments 
have indeed provided evidence for coherent vibrational 
dynamics[9| and ultrafast (typically ~ 100 fs) decay of 
electronic excitations. [lol Elf 

Of particular interest are ultrafast exciton dissociation 
processes at polymer heterojunctions, 0, 0, 0, 0, 0, 
comprising two materials with a AE off-set between the 
constituent valence and conduction bands (see Fig. 1). In 
molecular terms, this is the energy difference between the 
HOMO or LUMO levels of the two components. If AE 
is large compared to the exciton binding energy cb (for 
typical organic polymers, €3 ~ 0.5 eV), the photoexcited 
exciton state decays to an interfacial charge-separated 
state (or exciplex) . H U El El El Generally, such 
a situation is desirable for photovoltaic systems where 
AE provides the driving force for charge-separation upon 
photoexcitation. On the other hand, if eb > AE, the 
excitonic state is stable, and such situations are best 
suited for light-emitting diode (LED) applications. If 
€b — AE, charge separation and exciton regeneration 
processes tend to compete. 0, 0, 

Even though considerable progress has been achieved 
over recent years in the electronic structure characteri- 
zation of polymer semiconductor systems, pi. I5I Elf only 
tentative explanations have been given for the quantum- 
dynamical mechanism of the relevant charge transfer and 
decay processes. This applies, in particular, to the 



electronically nonadiabatic decay channels which are of 
key importance for the ultrafast phenomena mentioned 
above. Given the presence of multiple electronic state 
crossings in conjunction with electron-phonon vibronic 
couplings, nonadiabatic crossing regions and conical in- 
tersection topologies are expected to play a landmark 
role. Against this background, this Communication aims 
to provide a detailed quantum-dynamical picture of the 
nonadiabatic events determining exciton dissociation and 
formation of a charge transfer (exciplex) state at a poly- 
mer heterojunction. 

In the following, we focus upon the fate of the primary 
photo-excitation at a donor-acceptor heterojunction 
formed at the phase-boundary between poly[9,9- 
dioctylfluorene-co-N-(4,butylphcnyl)diphcnylamine] 
(TFB) and poly[9,9-dioctylnuorene-co-benzothiadiazole] 
(F8BT) when spin-cast from solution. El The 
optimized molecular geometry of the two polymer chains 
at the interface is shown in Fig. 1. LEDs fabricated 
using these materials exhibit considerable luminescence 
intensity. This is remarkable because the relative HOMO 
energy off-set between the two polymers, AE, is only 
slightly larger than e B • 0, II H [3 E3 The photoexcited 
exciton (XT) state is thus only marginally stable against 
dissociation into a charge-transfer (CT) exciplex state, 
and the majority of the primary exciton population 
is indeed found to undergo charge separation. El 0] 
Although the CT state is weakly emissive, appearing 
140±20 meV to the red of the exciton, the considerable 
luminescence is attributed to secondary excitons formed 
by endothermic back transfer of the hole from TFB+ to 
F8BT" to form the F8BT exciton with an activation 
energy of 100 ± 30 meV. El El El Nonadiabatic XT/CT 
state interactions clearly play a crucial role in this 
system. 

In order to understand how nuclear motion influ- 
ences the electronic transfer and decay mechanisms at 
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FIG. 1: (A) (Color online) Schematic diagram illustrating 
the electronic structure and photophysics of the TFB:F8BT 
donor- accept or heterojunction. The solid and dashed lines in- 
dicate the energy levels and band offset AE, respectively; the 
solid arrows indicate the electron occupancies. (B) Simpli- 
fied one-dimensional scheme of the photoexcitation and sub- 
sequent nonadiabatic potential crossing process. The black, 
blue and red curves indicate the ground, exciton and charge 
transfer states, respectively. The dashed lines indicate the 
adiabatic potentials at the avoided crossing. (C) Optimized 
molecular geometry, from density functional theory (DFT) 
calculations, of the parallel TFB and F8BT polymer chains 
at the heterojunction interface. 



the TFB:F8BT heterojunction, we carried out quantum- 
dynamical simulations based upon the two-state vibronic 
coupling model put forward by one of us in 

Refs.jlEi- 

A subpicosecond scale decay from the XT state is in- 
duced by vibronic interactions involving high-frequency 
(C=C stretch) and low-frequency (torsional) modes. As 
explained below, we identify an extended nonadiabatic 
coupling region in the vicinity of a multi-dimensional 
intersection space. While the electron-phonon coupling 
is largely carried by the high-frequency modes, the low- 
frequency modes are shown to play a key role in the de- 
cay dynamics. The main features of this dynamics are 



expected to carry over to a general class of conjugated 
polymer systems. 

The following analysis makes extensive use of a re- 
cently developed effective-mode description of the nona- 
diabatic dynamics at conical intersections, 0, E3- EiJ 
Efll EH ] in conjunction with efficient multiconfigurational 
quantum propagation techniques. E2 EE I3> EH Our 
previous analysis has shown that the cumulative ef- 
fects of many modes in a nonadiabatic coupling situ- 
ation can be represented by three effective modes that 
entirely determine the short-time dynamics. In the 
present study, we further connect this previous analysis 
to a systematic, Mori-chain type decomposition of the 
phonon bath[2(|E3 ( see a ls° the related development in 
Ref.[30|). As shown below, the explicit 24-mode phonon 




FIG. 2: (Color online) (A) Coupled adiabatic potential en- 
ergy surfaces (PESs) as a function of the branching plane co- 
ordinates (Xi,X2), and (B) zoom- in on the avoided-crossing 
region. The black, white, blue and red circles indicate the re- 
spective locations of the Franck-Condon (FC) geometry, the 
conical intersection (CI), and the minima of the exciton (XT) 
and charge transfer (CT) states, respectively. The blue and 
red lines indicate the XT and CT diabatic states, respectively. 
The dashed line indicates the XT-CT avoided-crossing seam 
line. 
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bath of the model under discussion can thus be replaced 
by a 9-mode effective bath over the time scale of interest. 

Our starting point is the following electron-phonon 
Hamiltonian for the coupled XT and CT states inter- 
acting with an A-mode phonon bath.[l4l Ejf 



N 



H 



= V A + J2 H * 



where 



H, 
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with pf = — d 2 /dx 2 and the electronic splitting Va = 
—Acr z (with cr z the Pauli matrix). Mass and frequency 
weighted coordinates were used, along with the conven- 
tion H = 1. As mentioned above, the phonon bath com- 
prises a high-frequency branch composed of C=C stretch 



modes and a low-frequency branch of torsional modes. 
An explicit representation in terms of iV=24 modes was 
introduced according to Ref . Q . 

Eq. Q corresponds to a general linear vibronic cou- 
pling form, in a diabatic representation,!^ and was pa- 
rameterized by semiempirical calculations as described in 
Refs.0, |]| 0- The Hamiltonian Eq. (Q) gives rise to an 
(A — 2)-dimensional conical intersection space, [2^] which 
is the key feature determining the nonadiabatic dynamics 
of the system. Indeed, the nonadiabatic decay induced 
by conical intersections is often ultrafast, i.e., occurring 
on a subpicosecond to picosecond time scale. 

For the purpose of the following discussion, a decompo- 
sition of the phonon bath in terms of collective, or effec- 
tive modes is introduced. To this end, we use an orthog- 
onal coordinate transformation, X = Ta;. [l7ill^ . ll^.l20| 
leading to the following decomposition in terms of effec- 
tive vs. residual modes, 



H = H c ff 



(2) 



J 



with the part H c h which contains three effective modes (X±, X2, A3) that entirely define the coupling to the 
electronic subsystem, 



H 



cfr 



Va 



6 o 



^KiXil 

i=l 



Xf)l+ dij(PiPj + XiXj) 1 

i,j=ij>i 




(3) 



Here, a topology-adapted representation^^ was chosen, where (Xi, X2) lift the degeneracy at the intersection (and 
thus span the branching plane [29lp. while A3 lies in the intersection space. The modes (Ai,A2) define a suitable 
reduced representation of the intersecting surfaces, as illustrated in Fig. 2. The parameters of Eq. Q relate to the 
original Hamiltonian Eq. (JTJ as described in Refs.jliLEil. 



The residual Hamiltonian H Tes contains the remaining 
(A — 3) modes, and their bilinear coupling to the effective 
modes and among each other, 

N O 

i=4 

N N , v 

»=1 j=4, ^ ' 

Importantly, H ICS is diagonal with respect to the elec- 
tronic subspace. 

The transformation leading to Eqs. HJ)-lQ} introduces 
a hierarchical structure in the phonon bath: While the 
effective modes couple directly to the electronic two-level 
system, the residual modes couple in turn to the effec- 
tive modes. The new choice of coordinates has important 



implications for possible approximations, since the effec- 
tive mode Hamiltonian J? e fl by itself preserves the first 
three moments of the overall Hamiltonian. 19] Hence, the 
short-time dynamics is entirely described by H c g. 

The concept of a mode hierarchy can be carried fur- 
ther by introducing additional transformations within 
the subspace of residual modes. In particular, H rcs can 
be transformed to a band-diagonal form, [2^ correspond- 
ing to a hierarchy of residual bath modes (see also Ref. [3(3 ] 
for a related analysis). An approximate nth-order, i.e., 
(3 + 3n)-mode Hamiltonian can be defined by truncating 
this hierarchy at a given order, 

n 

H^=H eS + J2HH (5) 

l=i 
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with the Ith order residual bath Hamiltonian 



3i+3 



i=3l+l 

31+3 i-1 , s 

+ E E diJPiPj+XiXAl (6) 

i=3l+lj=i-3 ^ ' 

The results shown in Fig. 3, to be discussed in detail 
below, illustrate the effect of truncating the hierarchy of 
residual bath modes. 

The scheme of Eqs. J5J-|J5J corresponds to a generalized 
Mori chain, which can be shown to conserve successive 
orders of the Hamiltonian moments. If truncated at 
a given order, the hierarchy can be formally closed by 
adding Markovian dissipation. [2|| I^J 

As can be inferred from the reduced representation of 
Fig. 2, in the (Ai,A 2 ) subspace, a conical intersection 
indeed exists in the present system, and an extended 
avoided crossing region is present in its vicinity. The 
figure further illustrates that the minima of both elec- 
tronic states are located on the same side of an avoided- 
crossing seam that departs from the conical intersec- 
tion. Following the terminology of Marcus theory, the 
system can thus be classified in terms of the "inverted 
regime". Q Due to the prevalent weak coupling in the 
avoided-crossing region, the zeroth-order picture of the 
dynamics is expected to be markedly diabatic. With a 
Franck-Condon (FC) initial condition (X\ = X 2 = 0), 
the wavepacket rapidly accesses the nonadiabatic cou- 
pling region, but does not directly encounter the conical 
intersection. 

Fig. 3 (trace "exact") shows results of quantum- 
dynamical calculations for the overall 24-mode sys- 
tem according to the Hamiltonian Eq. (1), using the 
multiconfiguration time-dependent Hartree (MCTDH) 
method. [23, 123L I24I [25} From the figure, one can infer 
that the XT state population has decayed to about 50% 
after 200 fs, with a corresponding increase in the CT state 
population. Subsequent to the initial decay, the popula- 
tions of the two states remain approximately equal, with 
an oscillatory behavior which features the characteristic 
period of the slow (torsional) modes (T tors i on ~ 300 fs) . 
Superimposed is a weak oscillatory structure due to the 
fast (C=C stretch) modes (T stretch ~ 20 fs). The dy- 
namics apparently remains in a coherent regime over the 
observation period, and does not reach an equilibrium 
state. 

The effective mode analysis as outlined above is ex- 
pected to shed further light on the mechanism of the 
nonadiabatic decay of the XT state. Fig. 3 illustrates 
the results of wavepacket propagation for the successive 
orders n = 0, . . . , 2 (i.e., 3 to 9 effective phonon modes) of 
the effective-mode hierarchy. These successive orders fea- 
ture in alternation the high-frequency modes (for n = 
and n = 2) and the low- frequency modes (for n = l).j2|| 



Several observations can be inferred from Fig. 3: (i) At 
the order n = (3 effective modes), i.e., for H = H e g, 
no decay of the exciton state is observed. While the 
effective-mode construction guarantees that the very ini- 
tial dynamics is correctly reproduced, 0] deviations from 
the exact dynamics occur from about 80 fs onwards, (ii) 
At the order n = 1 (6 effective modes), the XT — > CT de- 
cay over the first 400 fs is correctly described, (iii) At the 
order n — 2 (9 effective modes), the dynamics is correctly 
reproduced over the entire time scale of observation (1.5 
ps). 
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FIG. 3: Time-dependent population of the exciton (XT) di- 
abatic state, from MCTDH calculations for the overall 24- 
mode system (bold line, labeled "exact") and the successive 
effective-mode approximations if"' ' = iT e a (3 modes), if' 1 ' 
(6 modes), and ff' 2 - 1 (9 modes), see Eqs. (5)-(6). Note that 
the 9-mode result is very close to the exact 24-mode calcula- 
tion. 

Strikingly, the presence of the high-frequency modes 
alone, at the level of the n = description, does not 
account for the XT — > CT transition. A more detailed 
analysis shows that the dynamics is characterized by os- 
cillations between the adiabatic states, but the popula- 
tion essentially remains confined to the XT branch. The 
overall picture in the reduced 3-mode space is thus pre- 
dominantly diabatic. 

At the order n — 1 (6 effective modes), the low- 
frequency (torsional) branch of the phonon bath is added. 
As illustrated in Fig. 3, the initial XT state decay is 
now reproduced correctly. Apparently the low-frequency 
modes play a key role in the nonadiabatic dynamics, 
even though they do not directly couple to the elec- 
tronic subsystem in the effective-mode picture. A study 
of the energy flow among the effective modes reveals 
that the low frequency modes act so as to dissipate the 
energy contained in the 3-mode (Xi,X 2 ,X 3 ) subspace, 
by an intramolecular vibrational redistribution (IVR) 
process. |26j In addition, the non-adiabatic dynamics as 
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such is influenced by the low-frequency components. 

The next higher level, n — 2 (9 effective modes) , yields 
very good agreement with the 24-mode reference calcu- 
lation over the whole propagation period. The _H^ 2 ' ap- 
proximation thus provides a suitable surrogate Hamilto- 
nian on the time scale of observation. 

Finally, Fig. 4 shows the time-dependent [Xi,X2) po- 
sition expectation values for the XT vs. CT portions of 
the wavepacket, from the 9-mode calculation. A con- 
certed, oscillatory motion persists over the complete ob- 
servation period. The comparatively regular behavior 
is consistent with repeated passages through a weakly 
avoided crossing region, rather than a direct passage 
through a conical intersection. The presence of the con- 
ical intersection space is essential, though, in that it de- 
fines the relevant nonadiabatic coupling region. 

In summary, the present study emphasizes the role of 
coherent, quantum dynamical evolution, and its intrin- 
sically multi-dimensional nature, in characterizing the 
nonadiabatic decay of the photochemically accessed XT 
state. Note that the presence of a high-dimensional inter- 
section topology is not an exception, but rather the rule 
in multidimensional systems. |28l l3l| The effective-mode 
analysis proposed here leads to a reduced dimensional- 
ity description, which provides important insight both 
from a static and a dynamical viewpoint. The represen- 
tation of the potential in Fig. 2, in terms of the effective 
modes (Xi,X 2 ), yields a unique picture of the topology of 
the coupled electron-phonon system. From a dynamical 
point of view, the effective modes (X\, X 2 , X 3 ) determine 
the shortest, initial time scale, i.e., several tens of fem- 
toseconds. For the present system, it is however essential 
to (at least partially) include the residual (N — 3)-mode 
phonon bath. As demonstrated in Fig. 3, a 9-mode repre- 
sentation yields a very good approximation of the overall 
24-mode dynamics. 

One of the main conclusions of this study is that 
the presence of the low-frequency (torsional) modes is 
a key element in the observed decay dynamics. Even 
though the high-frequency (C=C stretch) modes are 
most strongly coupled to the electronic subsystem, they 
do not by themselves induce the observed XT state decay. 
The coupling to the low-frequency phonon branch gives 
rise to an energy redistribution (IVR) which is essential 
in mediating the nonadiabatic decay. The importance 
of the low-frequency branch has been previously recog- 
nized in the context of the static absorption/emission 
spectroscopy of similar systems. 0, H2 

Our observations are in qualitative agreement with ex- 
perimental results providing evidence for exciton regen- 
eration at the TFB:F8BT heterojunction.[l [r| Even 
though our calculations do not account for temperature 
effects, and the influence of dissipation is limited to the 
explicitly included phonon modes, we expect that the co- 
herent, non-equilibrium nature of the dynamics remains 
a dominant feature over the first few picoseconds. Fur- 




FIG. 4: (Color online) Time-dependent position expecta- 
tion values of the XT (blue) and CT (red) portions of the 
wave-packet as a function of the branching plane coordinates 
(Xi,X2), for the n = 9 calculation. The black, white, blue 
and red circles indicate the locations of the FC geometry, 
the conical intersection, and the XT vs. CT minima, respec- 
tively. The dashed line indicates the XT-CT avoided-crossing 
seam line. (A) Projection on the (Xi,X%) plane and (B) 3- 
dimensional (Xi,X2,t) plot. 

ther studies including temperature effects will allow for 
a detailed comparison with the results, e.g., of the non- 
Markovian master equation calculations reported in Ref. 
[0]- 
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